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Abstract 

The revised Enskog approximation for a fluid of hard spheres which lose energy upon collision is discussed for the case 
that the energy is lost from the normal component of the velocity at collision but is otherwise arbitrary. Granular fluids 
with a velocity-dependent coefficient of restitution are an important special case covered by this model. A normal solution 
to the Enskog equation is developed using the Chapman-Enskog expansion. The lowest order solution describes the general 
homogeneous cooling state and a generating function formalism is introduced for the determination of the distribution function. 
The first order solution, evaluated in the lowest Sonine approximation, provides estimates for the transport coefficients for the 
Navier-Stokes hydrodynamic description. All calculations are performed in an arbitrary number of dimensions. 
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I. INTRODUCTION 



Hydrodynamics is important because it provides a connection between microscopic models of particle interactions 
and experimentally observable behavior. Its experimental significance is due to the fact that hydrodynamics describes 
the evolution of fundamental quantities, including local mass, momentum and energy density which are important in 
a variety of applications from microfluidics to cosmology. At a microscopic level, these quantities are " slow" variables 
which evolve on a timescale which is well separated from faster, microscopic variables so that the effect of the latter 
can be adequately encapsulated in the various transport coefficients appearing in the hydrodynamic description. 
Recently, much interest has centered on the hydrodynamics of granular media which are characterized by a loss of 
energy during collisions between the grains 

0, Hi U- The simplest microscopic model of granular materials consists 
of hard-sphere grains which lose a fixed fraction e of the part of the kinetic energy associated with the longitudinal, 
or normal, component of the velocities at the moment of collision while the transverse components of the velocities 
of the colliding particles are unchanged. The transport properties of this model system, which will be referred to 
as a "simple granular gas" below, have been worked out in some detail [J, However, there is interest in more 
complex models in which the fractional energy loss e is itself a function of the normal kinetic energy as such models 
are apparently more realistic @)0j@- Furthermore, there are a wealth of phenomena, such as endo- and exo-thermic 
chemical reactions, in which kinetic energy gets converted to some other form and thus couples, e.g., chemical reactions 
and hydrodynamics. A good example is somnoluminescence[^| . [Toj where classical hydrodynamics is often used, with 
some success | llj. to try to understand the complex processes taking place under extreme conditions. Recently, the 
possible scattering laws for dissipative collisions consistent with conservation of momentum and angular momentum 
has been discussed and formally exact expressions for the balance equations of mass, momentum, energy and species 
have been formulated and the Enskog approximation discussed [12j. The purpose of the present work is to show that 
the Chapman-Enskog expansion ^3] can be applied to that kinetic theory for the case of a one-component fluid in 
D dimensions with an arbitrary model for the normal energy loss and reduced to a relatively simple form thereby 
providing convenient expressions for the transport coefficients covering this entire class of models. These results will 
therefore include as special cases the known results for elastic hard-spheres in 2 and 3 dimensions^) jU3 an( i f° r the 
simple granular gas in 3 dimensions^) as well as such interesting models as those with velocity-dependent coefficient 
of restitution0,[3,|3 for which the Enskog equation has not previously been solved. 

The emphasis on instantaneous interactions is due to their unique properties. Instantaneous hard-core interactions 
can be described in the Enskog approximation which is a finite density approximation known to give reasonable results 
for moderately dense fluids |l3l| and which can even describe hard-core solids^]. For most other interactions, only the 
Boltzmann description, a low density approximation, is available. The fact that the Enskog kinetic model is applicable 
to solids means that there is scope for application to dense granular systems, which might involve jamming, as well 
as to extreme conditions such as occur during somnoluminescence. Thus, despite its artificial nature, the hard-core 
model is an important tool in understanding the real world. 

In Section II of this paper, the possible scattering laws, exact balance laws and consequent Enskog equations 
are reviewed. The general form of the Chapman-Enskog expansion is also discussed. The zeroth-order Chapman- 
Enskog result is just the exact description of a spatially homogeneous system and is discussed in Section III. For an 
equilibrium system, this is the Maxwell distribution but when there is energy loss, the fluid cools and the resulting 
homogeneous but non-stationary state is known as the Homogeneous Cooling State (HCS). Section III formulates 
the description of the HCS in terms of a generating function formalism and gives all information needed to calculate 
the simplest corrections to the Maxwell distribution for arbitrary energy-loss models. In Section IV, the Chapman- 
Enskog expansion is extended to first order which is sufficient to get the Navier-Stokes transport properties. The 
general formalism is illustrated by application to the simple granular gas in D dimensions. The paper concludes with 
a discussion of the approximations made and comparison to more complete calculations. 

II. THEORETICAL BACKGROUND 
A. Kinetic theory with energy loss 

We consider a collection of particles having mass m and hard sphere diameter a which interact via instantaneous 
collisions. The position and velocity of the i-th particle will be denoted by ~q \ and v i respectively and the the 
combined phase variable ( q j, v %) will sometimes be denoted as cCj. The only scatte ring law allowing for energy loss 
that is still consistent with conservation of total momentum and angular momentum, 12j is that two particles having 
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relative velocity v 12 



v 1 



v 2 prior to collision must have relative velocity 



"' 7 I2 <j\2 j T u ■ </|2 Cr'u ■ 1] 1 2 I \/ (~V12 ■ §12) 2 - — 



(1) 



after the collision, where gx2 is a unit vector pointing from the center of the first atom to the center of the second 
atom. The energy loss is SE which we allow in general to be a function of the normal relative velocity 



SE = A ( u\ 2 • q) ■ 



(2) 



and the center of mass velocity V 12 = v 1 + u 2 is unchanged. It is easy to confirm that the change of energy upon 
collision is 
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The simple granular fluid model is based on the energy loss function 



A(x) 
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(4) 



where the constant a is the coefficient of restitution. More complex models involve a velocity dependent a while 
other choices of energy loss function, involving e.g. a thresholds for energy loss, would be appropriate for chemical 
reactions. In the following, it is convenient to introduce the momentum transfer operator b defined for any function 
of the velocities g uj;t) by 



bg{~v \,!?2\ t) = ffOw'i, ~v" 2 ;t) 



(5) 



In some applications, the energy loss may not occur for all collisions but rather might be a random occurrence. We 
will therefore also consider throughout a somewhat generalized model in which for any particular collision, the energy 
loss function is randomly chosen from a set of possible functions {A a (x)} with probability K a {q\2 ■ w 12) which may 
itself, as indicated by the notation, depend on the dynamic variable 512 • v xi- A simple case would be the a fixed 
fraction Ki = 1 — p of collisions are elastic with energy loss function A x (x) = while the remainder occur with 
probability K2 — V are inelastic with energy loss A2 (x) ^ 0. In any case, it is assumed that ^ a ^ a ( x ) = 1 f° r a ^ x - 
The momentum transfer operator for the type a collisions will be written as b a . 

The kinetic theory, Liouville equation and the Enskog approximation, for arbitrary energy loss function has been 
discussed in ref. |l2j| . The one body distribution function f(qi, Ui;t), giving the probability to find a particle at 
position ~q 1 with velocity v 1 at time t, satisfies, in the Enskog approximation, the kinetic equation 



d _> d 
ot a q 1 



f(x i; t) = J[fJ] 



(6) 



where the shorthand notation x% — is used. The collision operator is 



J If J] 



dx 2 T- (I2)x{~q i,~q 2; [n])f(~q 1, ~vv,t)f (~q 2, ~V2\t) 



where x(~qii~Q2] [ n ]) is the local equilibrium pair distribution function which is in general a functional of the local 
density. The binary collision operator T_ (12) is 



T_ (12) = 



Jp C"l>^2) (ba) K a (qu ■ ^12) - 1 



S (-l^ia • ~q 12) 5 (q 12 - a) ~vi 2 ■ qu 



(7) 



where (x) is the step function, (baj is the inverse of the momentum exchange operator and J a is the Jacobian 
of the inverse collision 



J a (^1,^2) = 



d 



(b a l?i,b a 1 V2) 

9(1?1,1?2) 



(8) 
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It will not be necessary to work much with this complicated operator as most calculations can make use of its simpler 
adjoint T+(12) defined for arbitrary functions of the phase variables A (xi, X2) and B (11, X2) by 



so that 



dxid,X2 A (x\, x 2 ) T_ (12) B (x\, X2) = — / dx\dx2 B (x\,X2) 2+ (12) A (x\, X2) 



^ K a (q 12 ■ ~vi2) b a - 1 



T+ (12) = 6 (- «\ 2 ■ ~<f 12) <5 (912 - o") (-V12 • gi 2 ) 



(9) 



(10) 



B. Hydrodynamic fields and balance equations 



The hydrodynamic fields are the local mass density p (V, t), the local velocity field u ( r , t) and the local temper- 
ature field T ( r , i). They are defined in terms of the distribution by 



p(~r* 7 t) = mn(r,t)=m / dv % f (r, v %; t) 



(11) 



p(r ,t) u(r ,t) = m J dvi v i/( r , v i;i) 

—n{~ 7 r 1 t)k B T{~r 1 t) = -m I dv 1 F x / ( r , v 1; t) 

where n(/r,t) is the local number density and D is the number of dimensions. In the third equation V \ = ~v\ 
~u(~r,t). Their time evolution follows from that of the distribution and is given bv|l2| 



d_ 

dt 



n + V • (Itn) = 



0_ 

dt 



p~u + V ■ (pltli) + V • P =0 



dt ^ ) Dnks 



P : Vu + V • g 



The pressure tensor is the sum of two contributions P = P K + P c with the kinetic part being 

^(7^)=™ / <n?i /j(T > ,l?i,t)V > i^i, 
and the collisional contribution being 

P c {~r,t) = -^r^y^ /" dxidx 2 gi2gi2 (gi2 • ^12) (g!2 - cr) 8 (-gi 2 ■ V12) 



xxC9i,"9 2;N)/(xi;i)/(x 2 ;t)i^a(?i2 -"^12) / dy 5 (~r - y~qi - (1 - y)~q 2 ) 

Jo 



■ \ - V \-i ■ q\2 - *<t» I ''12 ' '/I 2 ) \/ (1?12 ' 912) 2 - ($12 • I/12) 



Similarly, the heat flux has a kinetic contribution 



and a collisional contribution 
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r°(r,t) 



o- ^2 J dx x dx 2 qi2 ($12 • V i2 ) 5 (<?i 2 - cr) 9 (-g i2 ■ I/12) 



X X( 3 i, 92; [n])/ (a;i;t) / (x 2 ;t)K a (q 12 ■ v 12 ) / dx S(r - x q i - (1 - x) q 2 ) 

Jo 



X 2" + ' ( ~^ 12 ' ~~ s 9 n C^i2 ■ $12) Y (W*12 ' 9i2) 2 - ^-A Q (qi 2 • I/12) 



(12) 



(13) 



(14) 



(15) 
(16) 



Finally, because of the possibility of energy loss, the equation for the temperature includes a source term given by 
£C^>*) = 7^7 X! J dx i dx 2 ($12 • ^12) S (qi2 - a) 6 (-§i2 • ~v 12) (17) 

xK a (qu ■ ~vi2) A a (qi2 ■ ~vu) x (~qu~^2; N) / (xi]t) f (x 2 ;t) K a (x 12 ) 8 ( r* - if 1) . 

All of these expressions are exact, given the Enskog approximation, and show that the hydrodynamics of the system 
is completely specified once the one-body distribution is known. 

C. Chapman-Enskog expansion 



The Chapman-Enskog expansion is basically a gradient expansion of the kinetic equation assuming a particular 
form for the solution. Specifically, one attempts to construct a so-called normal solution in which all space and time 
dependence occurs through the hydrodynamic fields 



/(r*,u*;i) = /(v* |~r,V>t) 



(18) 



where the compact notation for the set of hydrodynamic fields ipt(r) — {n (r ,t) ,T (r ,t) (r , i)} has been 
introduced and the notation indicates that the distribution is a functional of the hydrodynamic fields at time t. This 
means that time derivatives will be evaluated as 



0_ 

di 



E 



d r 



dt 



T J(r,v-t). 



(19) 



To order the terms in the gradient expansion, we introduce a uniformity parameter e and replace V with eV and 
order terms in e. Since the space and time derivatives are related by the balance equations, we also introduce an 



expansion of the time derivative = d[ °^ 



as well as of the distribution itself 



f(~qi,l?i,t) = fa [v 1 \~q 1 ,ipt] + e/i [~vi \~q 1, ipt] + 



(20) 



Finally, in the Enskog approximation the collision operator is nonlocal and so must also be expanded (see Appendix 
as J [/, /] = J [/, /] + eJi [/,/] + .... Substituting these expansions into the Enskog equation and equating terms 
order by order in egives a set of equations for the distribution the first two of which are 



d q 



di Q) f ( Xl ;t) = J [/ ,/o] 
+ ) fo (xv,t) + df ] h (xi;t) = Jo [fo, fi] + Jo [fijo] + Ji [fo, fo] ■ 



(21) 



Since all time and space dependence of the distribution occurs via the hydrodynamic fields, the balance equations 
must also be expanded giving at zeroth order 



n(°) 
a t n 

q(0) -> 
o t nu 

g(0) T 



(22) 



Dnkr 



K0) 



and at first order 



d^p+v -(Hp) 







(23) 



'0, 



(i) 



u ■ V T 



o t pu 



Dnkf 



V • (puu) + V • P (V> = 

*p(0) 



V u + V • q 



(0) 



Dnk, 



(i) 



where, as noted, the fluxes and sources must also be expanded accordingly (see Appendix[5J). The logic of the normal 
solution is that these balance equations define the meaning of the time derivatives so that the time derivatives in 
ca. (|2f I are evaluated using ea. (|23|l and ea. (|f 9fl . Together with the expressions for the fluxes, eas. ljf 3H - 117I) suitably 
expanded, this gives a closed set of integro-differential equations for the distribution function. 
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III. CHAPMAN-ENSKOG AT ZEROTH ORDER: THE HOMOGENEOUS COOLING STATE 



A. Expansion of the zeroth-order distribution 

At zeroth order in the gradient expansion, the distribution /o {x\\ t) must be a local function of the hydrodynamic 
fields so eqs.(J2IJ and give 

( 2 m \ d 



yDn k J (, )^ f0iXl]t) = Mf ° jQ] (24) 
with 

C (0) ( V, t) = [ d^id^ 2 dq (q- ~v 12 ) 6 {-q ■ u\ 2 ) K a (aq- V 12 ) A a (a~q, ^12) (25) 

xxo (o-; n (~r)) fo {~r , «\;f) /o ( r*, "^2; ■ 

where xo ( <T ; n ) is the pair distribution function for a uniform equilibrium fluid of density n. Notice that this is not 
only the zeroth order component of the Chapman-Enskog expansion, but that it is also an exact (within the Enskog 
approximation) equation for the distribution of a spatially homogeneous system. If there is no energy loss, the solution 
will simply be the Maxwell distribution. When energy is lost in collisions, and in the absence of external forcings, the 
system cools and this is commonly known as the Homogeneous Cooling Solution (HCS). 

To solve for the HCS, we expand the velocity dependence about an equilibrium distribution by writing it as 

/o(*i) = hi (vi^t) £ c, {i, t ) S t (2^?) (26) 

where the Maxwellian distribution /m (vi; ^t) — nn~ D / 2 ( 2k J f l T ) exp (— 2 fc"r ' Ll i) ana - ^ ^ s important to note that 

this depends on the exact local fields ipt(r ). The functions {Si (x)}^ comprise a complete set of polynomials which 
are orthogonal under a Gaussian measure so that 



dlT hi (vi-^t) Si [n-fvi I Si I —vf 1 = AA, (27) 

where Ai is a normalization constant. In fact, these can be written in terms of the Sonine, or associated Laguerre, 
polynomials 

£ — ' 1 [a + m + 1) (k — m )\m\ 

rn—O v / \ / 



and satisfy 



so that in D dimensions 



' dx x a Lk (x) (x) exp (-as) = ^±A+ll 5mk (29) 

r (k + 1) 



and 



S k (x)=L k * (x) (30) 
T(iD + k) 

Ak r(iD)r(fe + i)- (31) 

Substituting eq.(|2SJl into the differential equation, eq.(|23J, multiplying through by L k 2 ( 2 fc™r ^i ) an ^ integrating 
gives 

(d^T^^) T ^f c ^^) + k(c k (A)-c k ^(yj t )) =^4, rs W<vWc,(W (32) 



G 



with 



h,rs (4>t) 



rn 



2k B T 

D-2 



vf T_ (12) 



2k B T 



-D 



hi {vv,ipt) 



(33) 



2k B T 



L, 



2k B T 



Notice that since £W ~ n 2 and Ik,rs ~ n i the coefScients Cfc can only depend on temperature so Ck (tpt) = Ck (T(t)). 
Furthermore, it must be the case that Co = 1 and c\ = in order to satisfy the definitions of the hydrodynamic fields. 
It is easy to show that I® a = so that the k = equation is trivial. Suppressing the dependence on ip t , the k = 1 
equation gives 



Dnk, 



(0) 



and it may be confirmed that this is consistent with ca. (|25|l . The first nontrivial approximation is to take c^. = for 
k > 2 and to use the k = 2 equation to get 



Dnk, 



_ I 1Q0 (7 1: 2o + h, 02) 



Dnk B T 



(0) 



+ 2c 2 



f'2 



^2,00 + (^2,20 + -^2,02) C2 



(34) 
(35) 



where the quadratic terms on the right are typically neglected as they are of similar structure to the neglected quartic 
terms. For a granular fluid, there is no quantity with the units of energy except for the temperature, so the coefficients 
of the expansion, which are dimensionless, are temperature-independent. For systems with additional energy scales, 
the coefficients must be determined by solving the resultant ordinary differential equations with appropriate boundary 
conditions. For example, if the energy loss were bounded, then at high temperatures it should be insignificant and 
one would expect limr^oo Cfc (T) = Sko to be the boundary condition. 



B. Generating function formalism 



The calculation of the integrals which define the coefficients on the left in ea . I|32|) can be formulated in terms of a 
generating function. Specifically, it is shown in Appendix IU1 that 



. r(ip)r(fc + i) (2k B T\ l/2 

h,rs (ipt) = -n -=rnr, 1 t\ — — T 
r (±d + k) v mcr / 



1 ,. ,. ,. d r d s d k 
x . .. . hm hm hm 



r\s\k\ zi^o z 2 ^o x^o dz[ dz% dx k 



J2 G * W*|A )-Go 



(36) 



where n* = ncr^.The generating functions are 



G a (ip t \K) 



-l ir - 1 /*S D (l-z 1 x)-$' 

/>oo 

x / du K* {y/u) exp 



1 — Z\X 





x exp 
x exp 



2 - x - z 2 - zi + xziz 2 

(2 — Z2 — Z\)x 1 



2 — x — Z2 — Z\ + xz\Zi 2 



1 - z 2 x 



2 — X — Z2 — Z\ + XZ\Z2 

1 (z 2 -z 1 )x 

22 — x — z 2 — zi + xz\z 2 



U — \/U\ u 



(37) 



and 



1 _i /n -D + l 1 

Go = — -7r 1 S B (1 — Zxx) 2 (2 — x — z 2 — z\ + xz\Zi) 2 . 



(38) 



7 



with Sd the area of the D dimensional unit hard sphere, 

Sd = T(Dj2) < 39 > 
which is, e.g., An in three dimensions and 2tt for D = 2. The scaled probabilities and energy loss functions are 



*"<•'•> M<\/^) (40) 



A* (x) = pA (xJ^) 
\ V mp J 



In general, K* (x) and A* (x) are functions of temperature and, hence, time but in order to keep the resulting 
expressions below from becoming too cumbersome, these arguments will be suppressed. The utility of this generat- 
ing function, which is admittedly complex, is that the limits and derivatives needed to evaluate ea. l|36[) are easily 
programmed using symbolic manipulation packages. 

To complete the description of the uniform fluid, I give the quantities necessary to calculate the lowest order 
correction. These are written conveniently as 

S D fk B T^ 1/2 



h ' rs = n * X 2D(D + 2)^ U^.' ILrs (41) 



roc 

Itrs = I*k E rs +J2 K a (y/Z) I^ s e-i U dU 

a J 

where Iq 22 + -^202 ~ ~& {D — 1) and all other elastic contributions are zero, and the inelastic kernals are 





T*I 

i l,00 — 




- 2) A* (y/u) 




J*I _ 


(A* 


(y/u) + 3 - u) A* (y/u) 


T*I 

i l,02 " 


r J l,20 — 


D + 
16 


2 (u 2 6u + 3) A* (Vu) 


T*I 

1 02,2 ' 


1 T*I _ 

r J 20,2 — 


16 \ 


' A* (01) (A* (u 2 - 
+16 (D- 



(42) 



The pressure can similarly be expressed in the generating function formalism, but here we just give for later use 
the expression for the pressure (see Appendix [Bj including the lowest order corrections to the Gaussian distribution 

Pi0) =l + n* X |£ (43) 



nk B T ^2D 



+nk B Tn*x^^=^2 1 vK * ("") 9 («> K (-«)) (l + ^2 (v 4 - 6v 2 + 3)) exp (-±v 
where the function g(v,A) is defined as 



dv. 



g (v, A) — sgn (v) \J v 2 — 2 A — v. 
C. The simple granular fluid in D dimensions 

For a simple granular fluid with constant coefficient of restitution one has A* (v) — (3A = (l — a2 ) \ v2 

and K* (x) — 1 giving 

G(i/j t \A) -» -it- 1/2 S d {1- z lX y^ D (1- zia:)' (2 - x - z 2 - Zl + XZl z 2 ) 1/2 (44) 



X 



1 I- 1 

— - (2 — z 2 — zi) x (l — a 2 ) + 2 - 2z 2 x + (z 2 — Z\) x (1 — a) 
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and 





T* — 
1 ,uu 




J 2,00 — 


^1,02 " 


' h, 20 = 


^2,02 " 


" ^2 : 20 = 







2 (D + 2) (1 - a 2 ) 
2 (1 - 2a 2 ) (1 - a 2 
3(L> + 2) 



(45) 



1 - a 



-8 (D - 1) + - (a - 1) (30a 3 + 30a 2 + 24Da + 105a + 137 - 8£>) 
8 

- (a + 1) (30a 3 - 30a 2 + 24Da + 105a - 56£> - 73) 



so that the lowest order correction to the Gaussian is 

^2,00 



C2 = 



— 2/1,00 — (^2,20 + 12,02) 

16(1 -2a 2 ) (1-a) 



24D + 9 - (41 - 8D) a + 30a 2 (1 - a) 

The cooling rate is 

e (0) = - (1 - a 2 ) n* X ^L (^) 1/2 nk B T 
and for a simple granular fluid, g (v, A* (— v)) = v (a — 1) gives 

p (a) =nfc B T jl + n* X |g(l + a)j 



16 



(>2 



(46) 



(47) 



(48) 



so it is seen that the second order terms do not contribute. 

This completes the discussion of the zeroth order Chapman-Enskog solution which is also the HCS. The simplest 
approximation is to take the distribution to be Maxwellian with the temperature obeying ea. (|22|l . The first nontrivial 
approximation is to keep C2 which is given by ea . (|34|l with coefficients that depend on the energy loss model. 



IV. CHAPMAN-ENSKOG AT FIRST ORDER: THE NAVIER-STOKES EQUATIONS 

The first order equation can be written as 

5 t (0) /i - Co [ft] = Ji [f , fo] - (d { t 1} + 1? ■ V) f (49) 

where £0 [fi] — (Jo [A 7 fo] + Jo [fo, fi}) is the linearized Boltzmann operator. The first order balance equations are 
used to eliminate the time derivative on the right. It is convenient to divide the first order heat source into two parts 

£1 = £0 [fx] + L [fo] (50) 

where the first term on the right is, as indicated, a linear operator acting on the first order distribution and the second 
is of first order in the gradients and depends solely on the zeroth order distribution. Furthermore, since it is a scalar, 
£1 [fo] must be proportional to the only scalar gradient, namely V • ~u so that we will write £1 [fo] — £,Y U [fo] V • ~u 
(see appendix IbI for more details). Then, the first order equation becomes 

d * 0) ^ + D^fZo [fi] (T A/o) -Co [/1] = Ji [fo, fo]-^&" [fo] (v • ^ (t A/„) -£ £ Bf (V; [fo]) 

(51) 
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with 



where the variable z 



2k B T 



V 2 



= U-7o + 



1 dpM 
nksT dn 



1 dp 



(0) 



nk B T dT 
2 



_d_ 

dzi 



fo 



9 f 
d 



DnkB 
m 



dT" 



JQ) 9 t mnV ? 
v gT Jo 2T 



fo V h 



2k B T 



V u Vi 



1 

D 



_d_ 

dz\ 
d 



fo 
fo 



It is shown in Appendix 1X1 that J\ [fo, fo] can be written as 



f 9 f 



Dnkf 



-f0 S; 



ainx (0)[ - 



(52) 



(53) 



where the detailed form of the operator Jr' is given in the Appendix. The right hand side of ea. (|51|l is therefore 
expressed entirely in terms of the gradients of the hydrodynamic fields so that the first order correction to the 
distribution must also be proportional to the gradients. Since the only vector available is V and the only tensors are 
the unit tensor and the symmetric traceless tensor ViVj — j^SijV 2 , the first order distribution must take the form 



fx = fo (xi) 



A("> (V i) Vud in + (v x ) V lt d t T + 

+y5V (aM) (?i) (VhVu - ^v 2 ) (a iUj + a 



Vi V • u 



(54) 



Then, both sides of the kinetic equation are expressed in terms of the gradients of the hydrodynamic fields and since 
those gradients can vary independently, their coefficients must vanish individually giving 



1 M/o]^/VL (Q) + I>7 



Dnk B T 



A {1) 



Co 



foA^^t 



(55) 



where Greek indices range over the four values n, T, Vu and du. In this equation, the capitalized index, I a , is a 
super-index corresponding to a set of Cartesian indices as illustrated by the definition 



*?„ (?) 



V h Vi,l, 



D 



D- 1 



1 

D 



(56) 



The linear functional K" encapsulates contributions coming from the action of the functional derivative on the non- 
local term VT as well as terms coming from [fi] and is given by 

2 



A^fo 



. d d 
~dn ^ ~d~T 

1 



DnksT 



to [fo. 



-Ao [fo] \foA { 



Vu) 



' DnksT 

The source terms on the right are, after some manipulation, given by 



dT Jo 



C n (?\f a 
C T (?\fi 



1 dpW 

n fo + 



nksT dn 



1 Sp(°) 



o = 



C Vu (?\fi 
C 9u (?\f 



nk B T dT 
2 



9 t 

dz 



9 f 
d~z f ° 



QT Jo 



= 



DnksT 
rn 



(^ V "[/o]-P (0)C ) (y/o 



mV 2 
2T 



2k B T 



9 f 

dz 



9 f 

dz 



(57) 



(58) 
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and 



fi?[/o,/o] 

[/o,/o] 

« Vn [/o,/o] 



^i 7 u [/o,/o] - 7 



1 ^ln n 2 x (o) 




(59) 



We conclude the discussion of the first order approximation with some general remarks concerning the solution of 
eas. (|55|) - (|59ll . First, the hydrodynamic fields are defined by ea. l|ll[l which can be written as 



dV 



(60) 



with the array of velocity moments ip (v^j 
zeroth order distribution satisfies 



A,i{r) = J il>i[y)f [T,V;t)dV 
so that it must be the case that all higher order contributions to the distribution give 



1, ^V 2 , V ) . However, from the definition ea. (|26|) . it is clear that the 



(61) 



0- j i>i(v)fj (?,V; t) 



dV 



(62) 



for all i and j. Since the gradients of the hydrodynamic fields are arbitrary, this means that in the case of the first 
order distribution, the coefficients of the gradients must be orthogonal to the first three velocity moments under the 

fo (V) 



measure 



or 



= 



fa (v) A™ (v) f 



V:t) dV. 



(63) 



Second, it is clear that ea. (|55|l is a linear equation in the coefficients [Vj so that the conditions for the 
existence of a solution follows the usual theory of linear operators. In particular, defining a Hilbert space with 



measure 



fo\Vj the Frcdholm alternative, which states that for linear operator L and source term B, the equation 

LV = B has a solution if and only if B is orthogonal to the null space of L. We expect that ip \ Vj is in the nu ll 

space of the operator defined by the right hand side of ea. H55|) . In fact, it is clear that multiplying by ipi (Vi) and 
integrating over velocities gives, on the right hand side, 



1 



DnksT 



/oj4 (Vn) 



d_ 

df 



(64) 



Now, the Cq term vanishes for ipi 
respectively. For the last choice, ipi 



1 and V due to the conservation of particle number and total momentum 



rn \r1 



V , it is only nonzero for a = Vu due to rotational symmetry (for other 



choices of a, 4>f a is a vector or traceless tensor). Thus the only non- vanishing element of this system of equations is 



that for a = Vm and ipi 



l -V 2 which becomes 



foA (Vu) 



dVi jV 2 Co 



foA (Vu) 



(65) 



and which is seen to vanish from the definition of [5], eq.{B23| and Co [<?]• Thus, ip is indeed in the null space of 
the linear operator and a necessary condition for the existence of a solution is that the right hand side is orthogonal 
to ip as well. That this is indeed the case is easily verified. 



11 



A. Approximate solution of the integral equations 

The integral equations summarized by ea.s H55|) - (|59[) will be solved by expanding the unknown functions in 
associated Laguerre polynomials as 



( 7 )=£ 



1 (a) r (- 2 2 2 + A °) 



2k B T 



V 



(66) 



where the coefficients are in general functions of the hydrodynamic fields, ai°^ — ai (ipt) although, for clarity, this 
dependence will be suppressed below. It is interesting to note that for a simple granular gas, Garzo and Dufty0 
wrote the first order distribution as in ea.l|53Jl but with f replaced by the Maxwellian Jm (ui;^>t). The use of / 
here is motivated by the fact that the source term in the first order equations, ea. l|55|) . is proportional to / so that 
it seems appropriate to use this in the definition of the first order correction but this is not necessary. In fact, one 
clearly has that 



fo(v)A^(v) = f (v)j2^ a) L 



(a) t { R 2 2l + X -) 



in 



2k B T 



V 



(67) 



2k B T 



V 



for new coefficients of which are linear combinations of the and in particular, if is the first non-vanishing 



coefficient in the sum, then ai^ = a),"' . In the development to follow, the usual approximation will be made wherein 
only the lowest non- vanishing coefficient is retained (the "lowest Sonine approximation") so that the two choices 
should be equivalent. However, in this approximation, the exact relation ai"^ = ai"^ is only preserved if terms of the 
form C2ai"' are systematically neglected. This therefore motivates the simplifying approximation made here whereby 
such terms are in fact neglected throughout. In the Conclusions, the present approximation is evaluated for the special 
case of a simple granular fluid in 3 dimensions by comparison to expressions given in ref. Q where all dependence on 
C2 is retained. 

In order to solve the integral equations, the expansions ea. (|66|l are substituted into ea. i|55|) and the a-th equation 

j). 2 ( 2fcflT ^i 2 ) an d tensorial indices contracted and V\ integrated. The left hand 

side of these equations is found to be simplified by the choices A n = At = 1, Av« = and Xg u = 2 which are made 
henceforth. The result after some simplification can be written as 



is multiplied by (j}f a ±j L 



'D 
v 2 



A„ 



k) (2k B T 



Dk B T{D/2)T{k + 1) V 



A Q r 



» 



dT 



where the contributions from the Boltzmann operator are 



— ka 
T 



(a) 
k-1 



I 



» 



= A5 



•«2 



(68) 



TO. 

1 kl 



dVi 4>t (Vi) if^ +Xa) 



2k B T 



L 



(^ + A Q ) 



2k B T 



V 



(69) 



and the last coefficient on the left is 



Ktf = J dV x cj> Ia (V,) Lp +A "» (^^) *5 
The source terms on the right are 



/„, [Vi)L 



A? 



dV 



2k B T 



(70) 



(71) 



and 



(^ + A Q ) 



2k B T 



vnnf [/ ,/o] 



(72) 
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A straightforward evaluation using the orthogonality and standard recursion relations of the associated Laguerre 
polynomials jig gives 



K kl 



T(\D + k + l) /2fcsT\ 
Y{\D)T(k + l) n \ m ) 



da'T^kl ( S an - V 8 a T-^ / ,, , 

on 01 J Dnk B 



Co 



(73) 



1 



2nk B T 



(0) 



foL 



(^) 



2k B T 



V 



— | ;; 

T 



and 



a; 



r(^ + fc + i) 
r(f)r(fc + i) 



2k B T (\ 



n 2k B T (\ 



(\D + k) 



1 rpdp 



!<>! 



+ 2fc + l) c fe - (fc + 1) c fc+ i - fcc fc _i + 



2 (tf*[/o]-P (0)c )(AKc*-i-<*) + r^O 



Set _ 9cfc + i 
dT 



Dk B T 



\ 



2nk B T I D-l (D+k)(D+k+2) , s 



/ 

(74) 

Finally, it is useful to note that the lowest order coefficients are related to the kinetic parts of the transport 
coefficients (see Appendix [5J). Specifically, the first order contribution to the pressure tensor takes the usual form 



PiP = -V (diUj + djUi - —Sij ( V • u 
where the shear viscosity, i], has a kinetic contribution given by 



V 



K 



-2nk B T 



k R T 



D 



7<% ( V • u 



d u 
l ) 



D - 1 

and the bulk viscosity 7 has no kinetic contribution. The the first order contribution to heat flux vector is 

~q (1) (~r , t) = -MVp - kVT 



(75) 



(76) 



(77) 



where n is the coefficient of thermal conductivity and \i is a new transport coefficient characterizing the way in which 
density gradients can cause heat flow due to differential cooling rates. It vanishes in the elastic limit. The kinetic 
parts of these transport coefficients are given by 



K 



,K 



nk B T 
nk B T 



k B T\ D + 2 



(78) 



fk B T\ D + 2 
1 m J 2 



These expressions are exact if f is replaced by a Gaussian in ea. l|54|l but if the first order correction is written in 

terms of /o then there are terms in C2 which would contribute (as discussed in the Appendix) in principle but which 

(7) 

would in any case be dropped here since they are of order c<2a So ■ The collisional contributions will be discussed below. 



B. Lowest order approximations 



The simplest nontrivial approximation is to keep only the lowest order nonzero coefficient in each expansion in 



ea. ()f)6[> . This means a\ 



(n) JT) (Vu) 



, u) and a^ 11 * 1 . Since the transport coefficients are more interesting than the distri- 



bution itself, we write these equations in terms of the kinetic parts of the transport coefficients giving 



1 T d\L K 
m dT 



£0 (D + 2) -T 

m 



9 In nxo 
dn 



nk B T 



k B T\ D + 2 



k B T D (D + 2) 



-C 2 

(79) 



2)i 

m 



T- 



k B T\ D + 2 



m 



1 m 4 
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l t D + 2 



T df a 



rVu JVu) _ n Vu 
1 22 u 2 — "2 



D + 2 (2k B T\ Br) 



r- 



<9T 



4fc s T 



(tf u [/oj - ;< 



(0). 



—10 



£>(£> + 2) -2 



L> - 1 



-ft 



da 



The Boltzmann integrals and the source terms are evaluated in a straightforward manner and the present evaluations 
were performed as described in Appendix[Dj making frequent use of symbolic manipulation The results can be written 
as 



Tl — tie 

J- r a t - a 



1 + E jf K (-«) e~^v (-v) S? s (v) + \vg (v, A* a (-„))) dv 



with the clastic contributions 



rnE jl E q ,. 



' li 



2D-1 , 



k B T\ 3/2 2(D-1) 



T VuE 

'22 



jduE 
J 00 



2 D-l c 

= ra a S D x 



Xn 2 a D ^S D 



k B T\ 1/2 (D-l) 
'k B T\ 5/2 AD 



rn 
'c b j 



and the inelastic kcrnals 



S11 (v) = ST(v) 



D + 8 



(v 2 - 1) 



1 



64 (D - 1) 
1 



64 (D - 1) 
— (« 2 - 1) 



16 (D - 1) 
(v 6 - 9v 4 + (8D + 49) v 2 - 37 - 8D) 

(w 4 - 6v 2 + 3) A* (u) 



The elastic contributions to the sources are 



ldn 2 X Da fk B T\ D + 5 
-a o~ 



2 dn 



-c 2 



2 D k B T 1 3 / dc 2 
rcr X& — " A I l + 2c 2 + — 



r?a D xSr> 



m T A 
D-7 



-n 2 a D X S D 



C2 

8D 

k B Tl ID -I 



m 2 V £> 



and the inelastic kernals are 
2$ («) 
T Q 9 o («) 



2 <9n 
k B 



16m 
1 

8£> 



(2 ((« 2 - 1) g (v, A; (-«)) + w 3 + 5«) A: (-«) - (9 - Av 2 + t, 4 ) 3 („, A; (-«))) 
((2A* («) + 3 - v 2 ) g (v, K (-„)) + «A: (-«) (« 2 - 1 - A* (-«))) 



1 / £>- 1 / fcgT 

2 V D I to 



(vA* a (-v)-g(v,A* a (-«))) 



(80) 



(81) 



(82) 



(83) 



(84) 
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and 



<?. M = i^^(^)i(« , -9« 1 + (49 + 8D)^-37-8D) a (»,A : (-»)) 

+^ (¥) a* <-> H 4 + 6 « 2 - 3 > (• <-» + "> 

UiAv) = -^—(-v 8 + Uv 6 + (-8D-88)v 4 + (126 + A8D)v 2 -2AD + 33)g(v 1 A* a (-v)) 

777 ZOD 

+ — 4i A : ("«) ((« 6 - n « 4 + 21 ^ 2 - 3) .9 (v, K (-„)) + « (t, 6 - 5« 4 + 9« 2 - 57)) 



m 128 
1 



128D 
1 



vA* a (-v) (A* a (-v) (-w 4 + lOv 2 - 15) + v 6 - lit; 4 + v 2 (61 + 8D) - 123 - 2AD) 



+ 128/5 5 (U ' A: ( 2A: ^ " ^ + 3 ) ~ ^ + 9v4 + ( 8C " 65 ) v2 - 8D + 53 ) 

° ' ' n ^ («A: (« 4 - lOt, 2 + 15) (t,, A: (-«)) (v* -6v 2 + 3)). 



The only non-vanishing coefficient of the temperature derivative is 

+ -La; (-„) ((« 4 - 6« 2 + 3) g (v, A* a (-«)) + (« 4 + 6t; 2 - 33) «) . 
The collisional contributions to the shear and bulk viscosity are, in this approximation, 



77° = -6rj K + 



D + 2 



7i 



7 = 7l -(nfc B T)ar-^=E H K v {3 - 6v 2 + v*) g (v, A* a (-v)) dv 

i2D\/ 2ir — Jo 



A' 



K+ = K K + ^ — 7l 



L>fc, 



2 TO 



— TO-xn (j 



2 D+l / K B 



fc B T' 



3/2 



1 S 



D 



m ) ATy/w D 



c 2 



1 + 1 E / ^ « 2e ^ 2 ( v2 - 3 ) 9 ^ A « (-«)) dw 



with 



71 = mtrn 



fc B T\ 2 S D n*x 



TO / D 2 y/lT 



1 1 

1 c 2 + - 

16 4 



(9 



3S7j 



-n X 



2D (D + 2) 

Finally, the first order corrections to the heat source are 

1/2 



^ ^ (-«) e"*» V 5 («, A: (-„)) (l + lea (« 4 - ICh; 2 + 15)) dv 

1 + ^ E J o K i (-«) ^" 2 (« 2 - 1) 9 («. A : (-»)) d « 



£i [/o] = (V • 1?) „^ X S D _*=_ ^ j( Kl (-«) A; (-„) w 2 e-*» 2 dt; 



+c 2 (V • 1?) n 2 a D X S D ^-=Y. jf K « ("«) ^ ("«) e " i " 2 ( 15 - + «*) * 



rfc 
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Equations Ip ^ -iP ty are the primary results of this paper. They give a prescription for the evaluation of the transport 
properties for an arbitrary model of energy dissipation at the Navier-Stokes level and in the usual, lowest Sonine 
approximation. In the next Subsection, these are illustrated by using them to give the transport properties of a 
simple granular fluid. 



C. Transport in simple granular fluids 

For the simple granular fluid, recall that A* (v) = (l — a 2 ) ^v 2 and g (v, A) — v (a — 1). Since there is no other 

;mperature as 

(91) 



(n) 


T -l/2 


(T) 


rp-3/2 


(Vti) 


rp-1/2 


a {9u) ~ 


rp-3/2 


Co ~ 


rp3/2 



giving 



& (2 3 + 2)- + iT l 



m 



dn ) 



nkgT 



k B T\ D + 2 



D(D + 2) k B T 

Q l + o Z~ C 2 



(92) 



£o (D + 2) 



Ml 



k k = mn 



k B T\ 2 D + 2 nT 1 2 fk B T\ 3 D{D + 2" 2 



(I + 202) 



3, D + 2 



i 22 



Vu \ „(Vu) _ qVu 



D + 2 /fc B T 



m 



+ / 9 "U* = mn 



fc R T 



J ~ 2fc B T 
2 



D(D + 2)-2 



D 



D - 1 



(9 If 



From the zeroth order solution, one has 



Co [/o] = - (1 - a 2 ) «*x 



, rr,\ 1/2 

k B T^ 



" 2y / 7T \m<7 2 

Equations l|8U|l are easily evaluated giving the Boltzmann integrals 



nk B T 



jn 

mi — 



tVu 

'22 



jdu 



k B T 



k B T\ 



3/2 



m 
1/2 



1 



1280F 



(a + 1) (3a (D + 8) - 1LD - 16) 
(a + 1) (30a 3 - 30a 2 + 105a + 24aL> - 56D - 73) 



(93) 



(94) 



n 2 a D - l S DX 



k B T 



5/2 



(1 + a) (3-3a + 2L>) 



and sources 
flj = 
= 



\^a D S D (^Pj 1 (1 + a) [6a (1 - a) - (3a 2 - 3a + 10 + 2D) c 2 ] 



m I 16 



(95) 



^-n 2 a D X S D — (a + I) 2 [2a - 1 + (a + 1) c 2 ] 
lb m 



n 



du 



6AD 
1 



rr<7 D x<SVj (1 + a) ( (5a - 1) (1 - a) (1 + a) - -c 2 (15a 3 - 3a 2 + 3 {AD + 15) a - (20L> + 1)) 



n 2 a D xSDk B T\ 



D-l 
D 



(a + 1) 



3a- 1 
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The low density (Boltzmann) transport coefficients in the elastic (a = 1) limit can be read off and are 

fk B T\ 1/2 ^ D(D + 2) 2 
K0 ^ kB {— 8a^S D D-l (96) 



Mo = 7o = 0. 

To facilitate comparison of the finite density transport coefficients to those of rcf . 4] , it is useful to introduce dimen- 
sionless quantities 

D-l, 3(£> + 8)(l-a)\ 

v; = x(l-^(l-a)(2D-3a-3) 

v* = —j-x (-4—) (30a 3 - 30a 2 + 105a + 24a£> - 56L> - 73) 



so that 



J_£(0) 



—;=bDnn 

\/7T 



m J \ ma 



rdu 



rVu _ 
J 22 — 



-=S D nn \ — \ f— 5 j „„ 



3 q *AsT\ 



1/2 



2 I ^7 



4-^/7? \m<r 

2D ,k B T (k B T\ 

^{D + 2) bDnn m [ma*) 



(98) 



and then 



n K = Kq 



D - 1 



D 



\v* T - 2(* 



l + 2c 2 



n*xS D (a + if (2a - 1 + (a + 1) c 2 ] 



8D(D + 2) 



fi K = 2k - [v* t - 2C\ 



, I K K* r ,D-l Idn^x 3(D-1) 

X \K Q H — — C 2 + - ^ ; JD 



(99) 



'/0 



D 

n -1 



^* - -C 



1 



2 <9n* 8D 2 (L> + 2) 
5c 



(1 + a) ^-2a (1 - a) + i (3a 2 - 3a + 10 + 2D) 2c^j \ 



AD{D + 2) 
i -l 



n*x (a + 1) (3a — 1) 



Vo 1 r * n^- x f *q \*^l^RiR±H (0)c 



where k k * = k k / kq and 
A* = {l + a) 



(5a - 1) (1 - a) (1 + a) - -c 2 (15a 3 - 3a 2 + 3 (AD + 15) a - (20L> + 1)) 

6 



(100) 
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The collisional parts of the transport coefficients are 

v o = n K S D f l + a\ _^ 

' 1 D(D + 2) A \ 2 J D + 2 1 K ' 

AS 2 D fl + a\ , 2 ( 1 



it(D + 2)D 2 V 2 / A V 16 

= 2D(D + 2) nX {—) fl 

« = ttt^ — t« X « + k (1 + a) — on X 1 c 2 

2D(D + 2) X \ 2 J 01 j 7r J D 2 GD + 2) 2 *V 16 

and the first order correction to the cooling rate is 

= 6>[/i]+£i[/o] (102) 

'3 , \ nk B T D + 2 



£o Lfi] = (v • if) ^ fl Tj (i - « 2 ) X 

a [/o] = ( V • t?) (|ti*bt) f§n* X (1 - a 2 ) 



Taking into account that the quantities ( = — 3 „ 2 flT £, c* = 2c 2 and cd = ^aj u are used in ref.jj, it is easy to 

verify that the present expressions agree for the special case of D = 3 with those of ref.Q up to terms of order c*aig 
except for the expression for 02 u ■ The expressions given for this quantity in ref.0 are incorrect and according to 
[l7j the correct expressions agree with those given here except for the terms — gc 2 (3 {AD + 15) a — (201? + 1)) in the 
expression for A* which give — |c 2 (81a — 61) for D = 3 whereas ref|lj| gives j^c* (159a — 19). 



V. CONCLUSIONS 



A normal solution to the Enskog approximation for a hard-sphere gas with energy loss has been determined using 
the Chapman-Enskog procedure to first order in the gradients, and the transport properties given to second order thus 
specifying the Navier-Stokes hydrodynamic description. The zeroth order distribution function, which describes the 
homogeneous cooling state, was expanded about equilibrium, ca. l|2()l) . and the equations for the coefficients given. The 
required collision integrals were expressed in terms of a generating function allowing for evaluation using symbolic 
mathematical packages and the explicit form of the required integrals needed to determine the first correction to 
the Gaussian approximation were given explicitly in the form of one-dimensional integrals. The expressions for the 
transport properties have similarly been reduced to simple quadratures for the standard lowest-Sonine approximation. 
The Navier-Stokes equations for such a system thus take the usual form, ea. (|12fl . with pressure tensor 

p i] = P (C % - V (diUj + d jUi - ^S i:j (V • 1?)^ - 7% ( V ■ 1?) (103) 

and heat-flux vector 

~q (~r, t) — — /i Vp — kVT (104) 

where [i represents a new transport coefficient not present when the collision conserve energy. Equations 179JI - 184J) 
determine the kinetic parts of the transport coefficients and ea.s H88f) and (|89|) determine their collisional parts. The 
pressure, p^°\ is given in ea. l|43() and the source term in the temperature equation, which accounts for the cooling, 
is given by eq.(@ZJ) and ea. (|9"0")) . Finally, as a simple application, the transport properties for a granular fluid in D 
dimensions were given and previous results for the special case D — 3 recovered. 

One question which has been left unanswered is whether we can judge the effect of neglecting some of the contribu- 
tions to the transport coefficients due to the non-Gaussian part of the zeroth order distribution. To answer this, we 
show in figure 1 the four transport coefficients for a simple granular fluid in three dimensions at a moderately high 
density of na 3 = 0.5 . The values are calculated based on the expressions of ref.jj] which include all contributions 
in c 2 , those given in the previous Section and the results of the Gaussian approximation, c 2 = 0. All three approxi- 
mations are in agreement for the shear and bulk viscosity, but the Gaussian approximation gives quantitatively poor 
results for thermal conductivity and the new transport coefficient. On the other hand, the expressions given above 
are in good agreement with the full expressions for the entire range of the coefficient of restitution thus giving some 
justification for the approximations used here. 
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FIG. 1: The four transport coefficients for a simple granular fluid in 3 dimensions at reduced density na 3 — 0.5. The lines are 
the results of ref.|^|, the circles are the results of this paper and the broken line is the Gaussian approximation. 
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APPENDIX A: EXPANSION OF THE COLLISION OPERATOR 

The collision operator is 

J[fJ] = -[ dx 2 T_ (12) X (91,92; N) / ( ~9i, ^V,t) f ( ? 2 , ^2; t) (Al) 

where x ( 9 1, ~Q2, [n]) is the local equilibrium pair distribution function which is in general a functional of the local 
density and we write binary collision operator T_ (12) as 



T_(T2) = 5{qi2-a)T_(12) 



T_(12) - 



J a (^1,^2) (bg) K a (qi2 ■ ^12) - 1 



6 (-^12 • 912) S (qi2 - a) ~vt2 • 912- 



Integrating over the argument of the delta-function gives 

J{f,f] = ~° D 1 J d ~v2<lq\2 t'_ (I2)x(~qi,~qi - vq\i\ N)/("? 1, ~vi;t) f(~q 1 - <rq\2,~v2-,t) . 

A gradient expansion of the nonlocal terms requires first an expansion of the one-body distribution 

/ On. - 0-912, ~«2;t) = / (~9 1,^2; *) - 0-912 • Vi/ (~q 1, l^a;*) + •■• 
For a normal solution, this becomes 



/On - 0-912, ~v2\t) = /(ti,^;*) - o9i 2 • ^2 ( v i^ On.) 



Hi ( 9 1) 



(A2) 



(A3) 



(A4) 



(A5) 
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and of the nonlocal dependence on the density of the pair distribution function 



(A6) 



Xo(<r;n(qi))+ / (fr (n(r)-n(q 1 )) 



8n(r) 



X (~Qi,~qi ~ ^$12; [n]) 



Xo (<r;n("?i)) + (Vin(lfi)J / dT* ("r - f 



Sn(r) 



which is accurate up to first order in the gradients. For a single-component system, it can be shown 1<§] that the 
second term reduces to a simple derivative giving 



x(~Qi,~Qi - 0-912; [n]) = xo (o-,n(~qi)) - ^ (<7<f 12 ■ V 1.71 ("if 1)) ^ X °^' — 



1(9*1) 



Using these results, we can write 



J [/, /] = Jo [/, f] + Ji [/, f] + 



with 



Jo [/,/] = -Xo / dx 2 T_ (12) f (~q 1 , If 1 ;t) f (~q 1 ,^u 2 ;t) 



(A7) 

(A8) 
(A9) 



which is, aside from the prefactor of xo (cj n {q i))i the Boltzmann collision operator. The second order term is 

5 



■M/./] = -^Xo^nC^i)) / dl? 2 dq 12 q 12 T'_ (12) 



dn(qi) 

which we write more compactly as 

•M/,/] = £ (vnfc (?i)) 

with 



( 9 1) 

<Tv 2 dqi2 (I12T- (12) / nfi,^?i;i) / ("if 1, V 2 ; t) 



/("?!, ¥ 2 ;tjjA10) 



5 



91nxo (cr;n("g i)) ^> 

"m ~ — /— > » J 1 [J; J 

on(qi) 



Ji[f,g} = / d~v 2 dq 12 q 12 T_ (12) f(~q 1 , ~vi\t)g(~q i,lj 2 ;t). 



(All) 



(A12) 



APPENDIX B: EXPANSION OF THE FLUXES AND SOURCES 

The expansion of the fluxes and sources is very similar to that of the collision operator described in the previous 
Appendix so only a few details will be given here. 



1. Pressure tensor 

The exact expression for the pressure tensor is P = P K ' + P c with 

*P K Cr,t\f) = m f d^ x f(J>,l? 1 ,t)ViV 1 , 



and the collisional contribution which can be written as 
m 

T 



P°(T,t\f) 



-v D ^2 J d!?id~v 2 dqqq (q- v*^) © (-<f • "^12) 



(Bl) 



(B2) 



xK a (q- v 12) / dy \ ( r + (1 - y)<rq, r - ycrq; [n])f(r + (1 - y) aq, v i;t)/(r - yaq, v 2 ;t) 
Jo 



x I - r l2 • q - san ( r ,2 • q) \ (1?i 2 ■ qf - ^A (q ■ ~v 12 ) 
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Clearly, the expansion of kinetic part is simply due to the expansion of the distribution function itself 



i=0 



with 



P*W (?,t)= P K (?,t\fi) 
The zeroth order contribution is based on fo which is homogeneous so 



i£ (r, t\f) = m-Sij I dV f ( F|Vt ) ^ = nk B T6. 



from the definition of the temperature. The first order contribution is 

h (V\i>t) ViVj 



P* (1) (r > ,t|/)=m / 



(B4) 



(B5) 



and comparison to the definition of the first order term eq.JSU shows that the only contribution is due to the shear 
term 



P l f 1) (~?Af) = m I dV f ( Xl ) 



''" k ~ £)^ feV ' U 



D 



dmk + 



\l rT~L A(9u) (?) ( VlVk ~ h SlkV2 ) { dlUk + dl 

diu k ~ ^Si k V ■ uj * + {SuSjk + S lk S jl ) J dV fo { Xl )A^ (v) 



ViV 3 {m) 



D - 1 

and using the expansion of A^- 9 ^ in associated Laguerre polynomials and their orthogonality relation gives 

1 



( r ,t\f) = mnSi 



(B7) 



> af u Cj 



2k B T 



-D/2 



uiuSpiT " -\l -jy—^ (diUj +d j u l - V • u 



exp (— x) L i 2 (x) L- 2 (x) x 2 dx 
1 



D(D + 2) 



4^ r(* + i) 



an i'2k B T 



(d - 2ci+i + c i+2 ) 



nk B T 



2a°» ( — )(! + <*) 



D 



D - 1 



diUj + djUi — j^Sij V ■ u 



and in the present approximation, we drop the term c%. 

The collisional part of the stress tensor will have terms arising from the expansion of the distribution as well as 
gradient terms arising from its nonlocality. The former gives a first order contribution of 

y c ^(^,t\f) = -j<r D Xo(<J;n(r))J2 I dlT 1 dlT 2 dqqq(q-ir i2 )Q(-q-ir i2 ) (B8) 

xK a ($12 • Vi 2 ) [fo Cr,^vi',t) fi { 1 r, 1 V2\t) + fi {~r,~vi 1 t) f (~r ,1? 2^)] 

x (^-~Vi2 ■ q\2 - sgn {~vi2 ■ $12) \j {~v 12 • qi2) 2 - ^A Q (q 12 ■ ^12)^ . 

while, using the results of the previous Appendix, the latter gives two terms 

tp c ( 12 )(y,t) = -^a D+1 X o(^n(y))J2 I dl? lC n? 2 dqqq(q-^ 12) e(-q-l? 1 2) (B9) 

xK a (q- ~v 12 ) -q - [( V/ ("r , "u fo (~r, w 2 ; t) - /o ("r , "u i;t) ( V/o ( r", "^2;*) 
x I -^12 • 9- sgn(~v 12 ■ q) J (tfu • <f) 2 - (<f ■ T? 12 ) J . 
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and 



ml D+i d Xo{^\n(~r)) 



dn ( r ) 

xK a (q- ~vi 2 ) f ("r ,l?i;t)/o ("r , V 2 ;t) 



Vin (7*)) • ^ / d~v idTu 2 dq qqq (q ■ ~v\ 2 ) 6 (— <?■ 7\ 2 ) (BIO) 



: ( - 7 , j ■ f/ - .sf//j ( 7* [2 ■ (/ ) i / (7*12 • qf - ^ A a (<f • ~v 12 ) 



so that C W = pCM+P^iP^ 13 ). However, it is seen that *P C ( 13 ' (~r ,*|/) = because of the integral is a 
vector but there are no zero-order vectors available from which to construct it. For similar reasons, P c ( 12 '> (~r, i)can 
be simplified to 

*P C ™{?,t) = -ja^Xoi^nir^id^Y, J dl? 1 d^ 2 dqqq(q-^ 12 )e(-q-^ 12 )K a (q-^ 12 ) (Bll) 



xqj 



k B T 



9_>_> _>_> _>_> d _> _> 

fo(r , «i;t)/o(r, v 2 ;t) - f (r , «i;t)F 2i — / (r, i; 2 ;i) 



x j <' .1.2 >/ - •"'.</" t r i> ■ (/) 5 ("?i 2 • qf - ^A a (q- ~v 12 ) 



where ^5^^. 



2. Heat flux vector 



The expansion of the kinetic part of the heat flux vector is treated analogous to that of the pressure tensor. It is 
given by 



qf (T, t\f) = ^ m I dir f (T, ir, t) vy\ 



1 



and the zeroth order contribution vanishes by rotational symmetry. The first order contribution is 
<zf (1) (^l/) - \m J di? fi WV 2 



(B12) 



(B13) 



Kn J d~v f (~r, V, t) AW (V) Vjdjn + A™ (V) VjdjT 
l —m j d~v f Q (~r, 77, t) [a™ (V) V A d t n + A^ ( V?) V 4 d t T 



2D 

where the vanishing contributions have been dropped. For both the density and temperature couplings, the important 
integral is 



l -m j erf f Q t) (v) 



2D 

1 1 (2k B T 

2 2D V m 



V 4 



(B14) 



V a^CjTT- / 2 / e~ x L7 (x) L~ (x) x^dx 



\s D ^mn (^pj 2 ]T ai a) Cj7 r- D / 2 e~* (if*" (x) - iff (x)) (l^ (x) - 2Ljff (x) + (*)) x^dx 



—-mn 
2D \ m 



T(A; + l)r(D/2) 



(§£> + 1 + 3fc) c fe - 2 (ID + 1 + ffc) c k+1 
+ (lD + l + fc) c fe+2 -fcc fe _i 
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So 



H K = [ nk B T 



knT\ D + 2 



(B15) 



K ^(„ t . sr (^)£±l) af ( 1 + (I) + 5)C2 _£±i C3 ) 



which gives the expressions in the text if the terms c 2 and C3 are neglected. 
The collisional part is 

~q C '(r,t) = -^° D J2 I d~vid-v 2 dqq(q- V12) 6 (-g • ~tf 12 ) (B16) 

a 

(g • ^12) / <fy X ( + (1 - y) <rQ, ~r - y<rq; [n]) f ( ~? + (1 - y) erg, f (~r — yaq, ~v 2 ;t) 

Jo 

x \( Vl+ V2 ) ' ^~^ 12 ' ^ ~ s 9 n (^ 12 • 5) y 12 ' ~~ ~A a ^12)^ • 

which gives at zeroth order 

~g c(0) (7*, i) = --^^Xo (o-;n(T*))^ y d~v 1 d~v 2 dq q{q ■ l?i 2 ) 6 (-g • l?i 2 ) 



(B17) 



xif a (g- V 12 ) /o ( r*,l?i;i)/o (7\l? 2 ;t) 
xi (?i+ ? 2 ) ■« 



r jo ,; - s.g// ( i ; L2 ■ 7) \/ ("« 12 ■ g) 2 - ^A a (g • "u 12 ) 



and at first order three contributions analogous to those of the pressure tensor 

-q C{11) (r,t) = -^<r D Xo(v,n(-r))J2 J d~v \dTu 2 dq q (q • V 12) 6 {-q • V 12 ) 

xK a (q- ~v 12 ) [f (7*, /1 (7*, r/ 2 ; t) + /1 (~r , /o (r*, V 2 ; t)] 



(B18) 



L2 • <? - s.gn- ("^12 • q) \j {~v\2 ■ qf - — A Q (g • ~v\ 2 ) 

m 



and 



>C(12) 



("?,*) 



m 
47 



(o-;n("r > ))y^ / dl? 1 dl? 2 dqq(q-~v 12 )B(-q-~v 



12) 



(B19) 



x^a (?• "u 12) ^g • (v/ (~r , /o (~r ,"u 2 ;t) - /o (~r i;i) (v/o (~r ,"u 2 ;t)) 



1 /- 



x- V"i + V 2 • g - u 12 • g - S5n(u 12 • g) \/ ( « 12 • g) A a (g • v 



and 



2 V 



1 ' j ~ 4^2 9n(V) 



m 



12 



^Vn("r)j y dVidV 2 dg gg (g • l?i 2 ) 6 (-g • V 



12) 



(B20) 



xif a (g • "u 12 ) fo (~r ,"v i;i) / ("r,^;*) 



+ V 2 ) ■ g ^- "U12 • g- sgn(~vi 2 ■ g) ^ (V12 • g) 2 - ^A a (g • ~u 12)^ . 

Now, the zeroth order heat flux vector vanishes by rotational symmetry of the homogeneous system. Only the 
vector parts of /1, i.e. those proportional to Vn and VT, can contribute to the first order contribution Ij ^ 11 ^. 
Similarly, the second term ~q c ^ 12 ^> can only depend on those gradients by obviously the contribution proportional to 
Vn vanishes. Finally, the third contribution vanishes as the integrand is odd in the total momentum. 
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3. Heat Source 

The heat source is 

= \^Y. I d^idlT 2 dq (q- V 12 )6(-g- V 12 ) (B21) 

a 

xK a (q- ~v 12 ) A a (?• V 12) X (~qi,~q 2 ; M) / (V, "u 1; i) / (~r - erg, ~u 2 ; t) 

so at zeroth order 

= \o D -\o(cr-,n{r))Y^ J dl! l dl?2dq (?• ^12) 6 (-9 • V 12 ) (B22) 

xif a (<? • 12) A a (§•■ 1? 12 ) /o ("r , /o (r*, 1? 2 ; t) • 

There are, as usual, three first order contributions but as in the main text, we separate these into £i(~r,t) = 
£o[/i]+£n("r,t) with 

Co [5] = ^^XofonC?)^ JdZiJrttdq (q- ^i 2 )S(-q- V 12 ) (B23) 

a 

xif a (g- V 12 ) A a (<?• V12) [/o (T',~vi 1 t)g(~r,lj'2;t)+g(~r,lj' 1 ;t)f (~r ,V 2 ;i)] , 

and 

£11 (r*,*) = ~a D '(drfi) / dvnTvidqqjiq- ~v i 2 ) & {-q ■ 1? i 2 ) (B24) 

a 

xif a (g- ~u 12) A a (<?• u\ 2 )x(~g\,~g 2 ; [n])/ ("r ,"u i;i) ^-/o (~r ,"u 2 ;t) 

and note that the coefficient of ax °(' T ^( r )) van j snes by spherical symmetry. The only nonzero contribution to incomes 

from tpi = ~u and the integral must be proportional to the unit tensor so we can write £n (~r, t) = £7" {~?,i)V -~u 
with 

H7 u (r,t) = -\°°Y. I dPJid!? 2 dq (q-l? 12 )e(-q-l? 12 ) (B25) 
xK a (q- v 12 )A a (q- v 12 )x(qi, 9 2; M) / ( r , ?; 1; t) / ( r , u 2 ;<). 



APPENDIX C: GENERATING FUNCTION FOR THE HCS 

The definition of the couplings is 



Irs,k = -n" 1 ^ 1 1 rfVid2 (2^?) T - ( 12 ) (^F) " (^1) ( C1 ) 

which is equivalent to 

Irs,k = n^Al 1 J d^^2 {^Pj /m(Vi)/m(1?2) (C2) 



• 1 - / in 2 \ ■' > - 



xL r 2 — » L s 2 — uf T+(12)L, 2 — vf 

\2k B T V V2fcsr V V2fc B T 1 
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The associated Laguerre polynomials are generated by 



1 d 
K (z) = -rlim 



n\ z^a dz n (l - z ) 



Q + l 



exp 



1 - z 



so that we have 
with 

Go 



m 



2k B T 



-1/2 



1 d r d s d k 

1 lim lim lim 



r\s\k\ zi^o z 2 ^0n o dz\ dz% dx k 



]Tg q (A o )-G 



1 



1 



-D 



(1- Zl )- (l-z 2 )~ (1-x)- 
x J d~v 1 d~v 2 dq cxp ^ ^ ^i - 1 ~^ 2^ " "^12) © (-§"■ ~« 12) exp ( _ j^T^i 



G a (A) 



1 



(l-z 2 ) T (l-s) 



11/"- 

d- !T 7r_£ ' / d ~"i d ~v 2 dqK* a {q-lj l2 ) 

Zo) * (1 - x) 2 J 



x cxp 



I-21 I-Z2 



v 2 (q- v 12) 6 (—q- v 12) 6 a cxp - 



l-x 



To evaluate the second quantity, note that 



v[ 2 = V 2 + V ■ ~v - V ■ q (l? ■ q + sgn (~v ■ q) \J (~v ■ qf - 2/3A (V • q)j + \ 2 - |a (~V ■ q) . 



A tedious calculation to complete the square in V gives 
1 1 1 



G a 



-D 



(l- Zl )~ (l-z 2 )- (l-x)- 
xK* (q- ~v 12 ) exp 
1 



dVdvdq(q- v)Q(—q - v) 



(2-z 2 -z 1 )x l A »/-> ^ 



x cxp 



1 

x exp [ - 



2 — x — z 2 — z\ + x Z \ z 2 2 
V 2 ) cxp 



1 x 
+ -, + 



1 — z\ 1 — z 2 1 — x 
(z 2 - z{) x 



1 — Z\X 



2 — x — z 2 — z\ + xz\z 2 



22 — x — z 2 — zi+ xz\z 2 



( ~v ■ qf + ~v ■ qsgn (~v ■ q) \J (~v ■ qf - 2A* (~u ■ q) ^ 



Performing the V integration and the D — 1 v integrations in directions perpendicular to q and finally the q 
leave the final result 



G a = --n- 1/2 S D (l-z lX ) 



-kr> 



1 — Z\X 



2 — x — z 2 — z\ + xz\z 2 



f-OO 

x du K* (y/u) exp 
Jo 



(2 - z 2 — z\) x 



x cxp 



x cxp 



1 

2 — x — z 2 — z\ + xz\z 2 2 
2 — z 2 x — Z\X 



2 \2 — x — z 2 — zi + xz\Z 2 
1 (z 2 -zi)x 



y/u-2A* a (V^y 



2 2 - x - z 2 — zi + xz x z 2 
We will also need the generating function for the case that no collision occurs 



Go = 



1 



1 



1 



-D 



(l-zi)- (l-z 2 )~ (1-z) 
x J d~v \d~v 2 dq exp 



v 2 - — ^ — v\ I q- ~v 120 (-<?• V12) exp 



1-zi 1 I-02 

Completing the square and performing the simple Gaussian integrals gives 



Go 



-7T 1/2 S D (1 - zxx) (2 - x - z 2 - Zl + xz lZ2 ) 2 . 
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APPENDIX D: EVALUATION OF INTEGRALS 



To illustrate the method used to evaluate the many integrals required in this work, consider the quantity 



►C(12) 



( r , t) defined in eq. (|B19|) and for convenience repeated here 



-$ c W(-?, t ) = -^a D+1 X0 (a;n(?))J2 I \dH 2 dq q (q ■ H 12 ) 9 (-? • H 12 ) 

xK a (q-lT 12 )^q- [(v/ (7^;t))/o(^V 2 ;0-/o(^,^i;^ (v/ (V, V 2 ; <)) 

(Vi + F 2 ) ■ q f- v*i2 • <f - ss" (^12 • 5) (^12 • qf - ^A a (q • v\ 2 )^ . 
As stated previously the only nonzero contribution comes through the temperature so we replace 

V/o - (VT) ±f . 
Keeping terms up to linear order in c 2 and defining the quantities 

m 

-V 2 



(Dl) 



2k B T 



(D2) 
(D3) 



V = 



2k B T 



2k B T 



( !?! + ~V 2 ) 



- ~v 2 ) 



we find that 



(^;/ (0) ^i;*)) / (0) ^2; t) - / (0) t) ^/<°> , t) 

(^/ (0) (r>, If i; *)) / (0) t) - / (0) Vi; t) ^/ (0) ^; t) 

UD 2 + ±D)-(D + 2)(V 2 + \v 2 



(D4) 



1 

T 



(-21 + z 2 ) 1 + c 2 



+ ( y 4 + ± W + + ( v ■ u 



i 

f 



2V -v \ l + c 2 



+ Zl c 2 (- ± (D + 2) + zi) - z 2 c 2 (-| (D + 2) + z 2 ) 
i(J) + 2) 2 -(fl + 4) (^ 2 + ^ 2 ) 



v 4 + ii/ 2 « 2 + ^> 4 + ( v • v 



/ M ( V, !m ( r", ^2; *) 



Substituting into the original expression and changing integration variables gives 

'2k B T^ 3/2 



qf™(?,t) = -^n 2 * D ^ X0 (a;n(r))(d 3 T) 



2k B T, 









1 + c 2 





xK a 

(V ■ q) [ -7 • q- sgn (If ■ q) \j (T ■ ,})' - 2 A;, j </ • T' ) 



7r D J]] / dVdvdq%qj (q • ~v) G (— g • ~v) 

a J 

|(D + 2) 2 -(D + 4) (F 2 + ±t; 2 ) 
f ( F 4 + \V 2 v 2 + ^ + (V ■ 1} 



(D5) 



Taking v to be the ^-direction for the V integral, this becomes 



C(12) 



n 2 a^ +i xo ((7;n(T))(a,-T) 



2fc s T 



3 /2 

/ dVdlTdqqSAQ-^fQi-q-^iDG) 



xK* a (q-l?)V x z l + c 2 



i(^ + 2) 2 -( J D + 4) (v 2 + ±0 
+ (y 4 + ±y V + ^ 4 + v> 2 ) 



-2V^ 



^-v* • q - sgn (~v ■ q) \J (~v ■ qf - 2A* (q- ~v)^j . 



2G 



The method now is to take q to be the cc-direction in the v integrals and to make the following substitutions: 

V 2 -» V 2 + Vl (D7) 

« 2 -> + 

and to integrate over V x , Vj_and v±. These integrals are performed by using symbolic manipulation to expand the 
kcrnal and to make make the replacement 



2tt 



to get 



,C(12) 



(r 



2 r // + l 



2" 



(D8) 



(D9) 



•K (v x ) Q + ^c 2 („£ - 2t£ - 9)) e->"2 - Sff n (« x ) 



m 
"17 



k B T\ 3/2 S 



x E u ' e ^ (l + ^ - 2v l - 9 )) (^-2A S ( Wx ) + v x ) dv x 

We can isolate the contribution due to A* (v x ) ^ by subtracting the A* (v x ) = term and using if* (w x ) = 1 
to get 



<1, '(r,t) = - T n 2 ^ +1 X o(a 4 T) 



3/2 



57 



1 7 1 

1 H c 2 + - 

16 4 



(D10) 



The only further complication is that some integrals involve kernals of the form h(V 2 ) (V -qj (V ■ 1?) . These 
can be handled by using the substitution 



(Dll) 
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